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Abstract 

The Helicopter Icing Consortium (HIC)f conducted one of the first U.S. 
tests of a heavily instrumented model in the controlled environment of a re- 
frigerated tunnel (lj. In the Icing Research Tunnel (IRT) at NASA Lewir 
Research Center (NASA LERC), ice was accreted on the main rotor blade 
of the BMTR-1 Sikrosky model helicopter under a variety of environmental 
conditions, such that liquid water content (LWC) and volume mean droplet 
diameter (VMD) ranges reflected the Federal Aviation Agency and Depart- 
ment of Defence icing condition envelopes (1). This report gives the corre- 
lated results of the data provided by NASA LERC. The method of statistical 
analysis is discussed. Lift, thrust, and torque coefficients are presented as a 
function of icing time, as correlated with changes in ambient temperature, 
LWC, and VMD. The physical significance of these forces is discussed. 

Introduction 

A majority of the U.S. military and civil helicopters in use today do 
not have clearance to operate in forecast idng conditions, or have a strictly 
limited clearance, due to the extreme sensitivity rotors display under ice 
accretion conditions. Currently the industry is lacking reliably adequate ro- 
tor deicing systems, and present all-weather certification of U.S. military 
and civil helicopters relies heavily upon flight testing in natural icing con- 
ditions, which is an expensive, potentially hazardous, and time-consuming 
process. Therefore U.S. helicopter manufacturers, in order to increase their 
all- weather capabilities, have expressed a desire for alternate methods of data 

This research is supported by NASA Lewis Research Center Grant 
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Graduate Reseach Assistant. 

f HIC is comprised of NASA Lewis Research Center, Texas A&M Uni- 
versity, Bell Helicopter Textron, Boeing Helicopter, McDonnell Douglas He- 
licopters, and Sikorsky Aircraft. 



acquisition to be used in obtaining FAA certification. Prior to the OH-58 
Tail Rotor Rig test at NASA LERC, only the French had conducted any 
icing wind tunnel testing of a scaled rotor (3). Their data compared favor- 
ably with full-scale natural icing flight test data, and indicated that icing 
tunnel testing would prove to be _a viable method of collecting iced rotor 
performance data (4, Appendix B). The promising results obtained in the 
ONERA Wind Tunnel by Guffond (3) was instrumental in the formation of 
the HIC, the OH-58, and the Sikrosky BMTR-1 testing. 

Data Acquistion 

Experimental data was obtained through a dynamic data acquisition 
system, a safety of flight system, and a number of other measurement and 
cataloguing techniques. Spray times, temperature and shedding event times 
were recorded, and each run was video-taped to provide a time-history of 
the ice accretion and shedding. Appendix C contains photographs of the ice 
accreted along the leading edge of the rotor blades during testing and lists 
the conditions under which the ice shapes were formed. 

Data Analysis 

Ratio of coefficients of torque, lift, and thrust to solidity, i.e., Cq/<t, 
Ci/c 7, and C-rja , respectively, were plotted versus icing time for correla- 
tion parameters of temperature, LWC, VMD. Early analyses of the data (5) 
indicated that shaft angle, 6, had little effect on torque, lift, and thrust re- 
sponse to icing conditions. Therefore two sets of data were correlated; one 
set includes 6, the other does not assume 6 to be a significant variable. 

Data files were received in Lotus format and plotted using the graphics 
package Grapher, readily available in the Department of Aerospace Engi- 
neering at Texas A&M University. The data was corrected for icing start 
time before plotting. In the original files, data was obtained before the ini- 
tiation of icing conditions, therefore when data was correlated, points were 
taken beginning at the icing start time. The last twenty to thirty seconds of 
data taken indicate the response of the rotor to a decrease in ice accretion 
as the spray rig was shut off, referred to as the rotor’s ’’post cloud behavior 
(2).” The standard deviation of the statistical variables LWC, VMD, and 
temperature was kept to a minimum when choosing data sets for analysis; 
however, no other restrictions were placed on the chosen sets. It was found 
that including 6 as a statistical variable greatly reduced the variety of icing 
runs which could be grouped for analysis. Effort was made to use as many 
different sets of icing data as possible in order to obtain a representational 
subset of the complete data taken. 


Physical Significance 

Helicopter rotor performance is highly sensitive to ice accretion; lift, 
thrust, and torque will vary with the type of ice accreted, the variation of ice 
growth along the span of the blades, and the amount and location of shed 
ice. Thus it is important to understand the characteristics of ice accretion 
and the physical properties of the accreted ice in order to begin to describe 
rotor performance variations. 

Ice will grow nonuniformly along the span of a blade as a function of 
local velocity variation and local aerodynamic heating effects. The smooth 
rime ice, which in some experiments has been shown to actually increase 
the performance of a blade under certain conditions, increases in thickness 
as a function of velocity. However, aerodynamic heating effects increase in 
significance as a function of velocity also, and can raise the surface temper- 
ature such that the rime ice shape is transformed into the more detremental 
double-horn shape of glaze ice (2). 

Rime ice is associated with smaller droplets, lower LWC values, and 
temperatures below 14 degrees Fahrenheit. These smaller droplets are more 
likely to be deflected around the airfoil, resulting in a small impingement 
zone near the leading edge. Glaze ice forms at temperatures between 14 
degrees Fahrenheit and freezing and has a density nearer to that of water 
than rime ice. As droplet size and density increases, the impingement limit 
moves farther back on the airfoil of the blade, causing a greater area to be 
influenced by the accretion process. 

The shedding process plays an important role in the outcome of an 
icing encounter. If ice is shed from a propeller late in the accretion pro- 
cess, the shed ice may cause damage to the fuseage or may be injested by 
intakes. Heavily accreted blades which shed nonuniformly may experience 
unbalance, which results in high vibratory loads on the rotor. The goal of 
deicing systems, then, is to force shedding early in the accretion process, be- 
fore significant decreases in lift and thrust occur, and before the mass (thus 
energy) of the shed ice becomes a hazard to the aircraft. 

Torque 

In all cases tested, torque rose as a function of icing time. For both sets 
of data, Cq/<t increases for increasing VMD, as shown clearly in Figure 1. 
Figure 2 shows a no-shedding case at VMD=15/xm. Rotor torque decreases 
slightly for minor shedding, and exhibits more pronounced changes when 
large areas of accreted ice are _shed. In almost cases, some residual ice will 
be seen after shedding. 


Increase in LWC increases Cq/<t as expected, as seen in Figures 3 and 
4. Note the overlapping which occurs in Figure 3 that does not in Figure 
4, suggesting that LWC plays less of a role at higher temperatures on the 
affects of torque. The effects of temperature at two different shaft angles are 
shown in Figures 5 and 6. As temperature increases, the torque needed by 
the rotor to maintain a contant ro t ational velocity decreases, indicating that 
some melting and/or shedding has occured on the blades. 

Figure 7 clearly indicates the effects of negative and positive 9 on torque. 
Both show 7 steady increasing functions with icing time, but the shift in value 
of Cq/( 7 is approximately 0.28 throughout the icing time, where Cq/<t{9 
=-6.5) > Cq/(t{9 =4-2.4). 

Lift 

Lift w'as shown to decrease with icing time for all parameters analyzed. It 
is not smoothly decreasing, but an overall oscillatory decrease. At the start of 
an icing encounter, lift will decrease rapidy due to the roughness induced by 
the initial ice accreted (2). Figures 8 and 9 show CL/cr as a function of VMD; 
there is much mixing of the data, indicating that lift is not strongly dependent 
on volume mean droplet diameter, and more dependent on other parameters 
such as liquid water content, which may have a greater influence on the 
weight of the accreted ice in normal icing encounters (excluding freezing 
rain). For increasing LWC, lift exhibits a more defined decrease, as seen in 
Figures 10 and 11, however even this decrease of LWC is only on the order of 
3 to 5% from LWC=0.35 gm/m 3 to LWC=0.75 gm/m 3 . In Figure 12, 
clearly increases over a temperature range of -12.3 °F to 4-27.4 °F. For the 
same change in 9 as in the case for torque, lift is not so clearly influenced; 
i.e., data overlap continually as Cl/< 7 decreases over icing time as seen in 
Figure 13. 

Thrust 

Ct/v exhibits the same trends as as expected, as lift and thrust 

are identical in nature for this case. Had forward motion of the aircraft been 
simulated then the values of lift and thrust would not have coincided, as a 
component of the total thrust affected by the ice accretion would have been 
delegated to forward motion only. Numerical differences between the two 
coefficients tend to enter in only the third or fourth decimal place. Figures 
14 through 17 show Ct/cf as a function of icing time for the parameters of 
VMD, LWC, and 9. See Appendix A for statistical data concerning plots of 
Ct/cf and Cl/ft. 


Concluding Remarks 

Examination of the correlated data has provided encouraging results. 
The data show the effects of liquid water content, volume mean droplet 
diameter, temperature, and shaft angle. Qualitative trends correctly follow 
those known to exist in flight. The model exhibits the correct physical trends 
for an actual icing encounter, including reactions to shedding and post-cloud 
behavior. Multiple runs with the s ame parameter values, though not plotted 
formally here, have shown good repeatability. Torque rises as a function 
of icing time for all parameters tested, while lift and thrust decrease. It 
appears that this data, and the method from which it was provided, will 
prove helpful in future research and development of model testing in the 
area of data acquisition for all- weat her flight certification. 
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Figure 4. Torque 
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Figure 9. Lift as a Function of Icing time 
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Figure 10. Lift as a Function of Icing time 
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Figure 12. Lift as a Function of Icing time 
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Figure 14. Thrust as a Function of Icing Time 
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Figure 15. Thrust as a Function of Icing Time 
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Figure 17. Thrust as a Function of Icing Time 
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Appendix A 
Scatter Plot Statistics 

Cl! a C t !<j 


AVG. 

ST. DEV. 

0.06375 

1.36E-3 

0.06394 

1.424E-3 

0.06045 

1.552E-3 

0.06190 

2.007E-3 

0.06238 

1.313E-3 

0.06267 

1.064E-3 

0.06178 

1.638E-3 

0.06135 

2.025E-3 

0.06182 

1.674E-3 

0.05894 

1.800E-3 

0.06116 

1.843E-3 


AVG. 

ST. DEV 

0.0631 

1.299E-3 

0.0636 

1.109E-3 

0.0604 

1.751E-3 

0.0619 

2.158E-3 

0.0624 

1.33E-3 

0.0628 

1.077E-3 

0.0617 

1.671E-3 

0.0615 

2.238E-3 

0.0618 

1.897E-3 

0.0590 

2.103E-3 

0.0611 

1.987E-3 


The above statistics corespond to the data in those plots for which a 
large spread in values appeared, making trends more difficult to read directly. 
Only those runs and figures noted here are represented above: Runs 70, 72, 
and 75 correspond to Figures 3, 10, and 15; runs 85, 90, and 91 correspond 
to Figures 1, 8, and 14; runs 99, 100, and 101 correspond to Figures 1, 8, 
and 14; run 87 and 89 correspond to Figures 7, 13, and 17. 
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ABSTRACT 

A consortium composed of NASA, Texas A&M 
University, Bell Helicopter Textron, Boeing Hel- 
icopters, McDonnell Douglas Helicopters, and 
Sikorsky Aircraft conducted an experimental pro- 
gram to Investigate the characteristics of a 
model rotor under Icing conditions. This pro- 
ject resulted In the first U.S. test of a heavily 
Instrumented model rotor conducted In the con- 
trolled environment of a refrigerated wind tun- 
nel, the NASA Lewis Icing Research Tunnel. The 
tunnel entry used a Powered Force Model with 
a 1.83-m (6 .00-f t)-di ameter main rotor, with 
0.124-m (4.9-ln. }-chord main rotor blades fabri- 
cated specially for this experiment. Test con- 
ditions Included a range of llauld water content 
and median volume droplet diameters that fell 
within the FAA and DOD Icing envelopes. The 
test data show the effects of icing on rotor 
lift, rotor torque, blade loads, and vibration. 
Ice shapes and Ice dimensions were taken, and 
molds were made of three Ice shapes. High-speed 
movies were taken to document Ice shedding. The 
results have been compared with analytical accre- 
tion predictions. 

NOTATION 

A c accumulation parameter 

B total pressure, kN/m 2 

C[_/o lift coeff Iclent/solldi ty 

C m model chord, m 

Cp full scale chord, m 

Cq/o torque coefficient/solidity 

Cj/a thrust coefficient/solidity 

c chord, m 

d droplet diameter, pm 

K 0 modified Inertia parameter 

LWC liquid water content, g/m 2 

R leading edge radius, m 

T static air temperature, Q C 

V free stream velocity, m/'se: 


a solidity (0.173) 

p advance ratio 

p ice density, g/m 2 

t Icing time, sec 

QR rotor tip speed 


INTRODUCTION 

The majority of U.S. military and civil 
helicopters have a restricted or no clearance 
for operating In forecast Icing conditions. 

This is due In part to the lack of adequate 
rotor deicing systems, extreme rotor sensitivity 
to ice accretion, and the expense of system 
qualification for the full Icing envelope. 
Procedures have been established to use Icing 
tankers and/or ground Icing facilities, but 
these approaches have not adequately simulated 
natural Icing conditions. Because of the trend 
toward design of all-weather rotorcraft, it Is 
necessary to develop and validate the experi- 
mental techniques to understand the effects of 
Ice on rotor performance and to determine the 
design requirements for rotor Ice protection. 
Analytical methods are now being developed to 
predict rotor performance In Icing, including 
Ice accretion and Ice shedding. The use of a 
model rotor to provide a less expensive and more 
repeatable source of test data wlli aid In the 
validation of these analytical methods and sup- 
port the overall goal toward more expedient 
certification procedures (Figure 1). The 
French undertook a model Icing program where a 
1/4-scale model rotor was tested In the ONERA 
SIMA wind tunnel at Modane, France (Refer- 
ence 1). This testing showed that a model pro- 
gram had promise, but the data acquired were 
too limited to provide a correlation base. How- 
ever, this testing prompted the Interest of 
rotorcraft icing researchers In the United 
States, A rotor Icing consortium evolved to 
Identify and pursue an approach to model rotor 
testing In an Icing environment. NASA Lewis 
Research Center undertook a program with univer- 
sity and Industry participants to demonstrate 
the usefulness of the model rotor test technique 
as an approach for obtaining meaningful icing 
data for rotating systems. 
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FIGURF 1 . - CERTIFICATION now DIAGRAM TOR ROlOnCHAfT ICING. 


TFCHNICAL APPROACH 

The Rotor Icing Consortium Is composed of 
members from the four major helicopter companies: 
Bell Helicopter Textron, Boeing Helicopters* 
McDonnell Douglas Helicopters, and Sikorsky 
Aircraft Division of United Technologies, plus 
representatives from NASA Lewis Research Center 
(LeRC), U.S. Army Aviation Research and Technol- 
ogy Activl ty-AVSCOM, and Texas A&M University. 

The consortium worked during the early part of 
the program to identify the necessary tasks and 
the means to achieve them. Each of the members 
was given specific responsibilities, as outlined 
i n Table 1 . 

A two model approach was selected as the 
most effective means to accomplish the program 
goals. A lightly instrumented OH-58 tail rotor 
(figure 2) that had been modified to approxi- 
mately operate as a main rotor was chosen as the 
initial test article. This would be followed by 
the more sophisticated Sikorsky Basic Model Test 
Rig ( BMTR) Powered Force Model (PFM) (Figure 3) 
to provide detailed rotor performance measure- 
ments. The OH-58 IRT entry would establish 
operational techniques for the more complex 
model, define tunnel test capabilities for a 
rotorcraft entry, and gather preliminary data. 
The OH-58 model test was completed in October 
1988 and the results presented at the 1989 Amer- 
ican Helicopter Society Annual Forum (Refer- 
ence 2). The conclusion of this Initial test 
indicated the tunnel was capable of being used 
for model rotor icing. The test techniques 
developed for the OH-58 entry were refined to 
provide a greater probability of success for 
the PFM entry. 

The Initial PFM Icing Research Tunnel (IRT) 
entry was scheduled for an 8 week period during 
July to September, 1989. During this time the 
model was Installed and icing testing began. 
However, blade tracking problems and the high 
humidity, warm amblent'condi t ions precluded 
testing at very cold temperatures. A total of 
44 icing encounters occurred in the 14 days of 
testing. 


At the end of this Initial test period, the 
model was removed and a second entry scheduled 
for November 1989. Between tests, the model was 
serviced, data reviewed, high speed motion pic- 
ture capability added, and methods for Improving 
blade track were prepared. This second entry 
resulted In 41 Icing encounters In 8 test days. 


MODEL AND FACILITY DESCRIPTIONS 
Icing Research Tunnel 

The IRT Is a closed-loop refrigerated wind 
tunnel. A 4100 hp fan provides airspeeds up to 
134 m/sec (300 mph>. The 21 000 ton capacity 
refrigeration heat exchanger can vary the total 
temperature from -1.1 to -42 'C (Figure 4). The 
spray nozzles provide droplet sizes from approx- 
imately 10 to 40 pm median volume diameter (MVD) 
with liquid water contents (LHC) ranging from 
0.2 to 3.0 g/r»w. The tunnel Is 1.8 m (6 ft) 
high and 2.7 m (9 ft) wide. 

Armor plates were attached to the walls of 
the tunnel test section, during the model rotor 
test, to protect personnel in the control room. 
Video systems were Installed to monitor the test 
area and local blade positions, and to provide 
tracking information. For several runs 0.635 mm 
(0.025 in.) aluminum sheets were placed on energy 
absorbing material and attached to the armor 
plating In line with the rotor disk plane. The 
shed Ice made permanent dents in the sheets that 
will be calibrated to calculate the Impact energy 
of the Ice that struck them. A high speed 16 mm 
camera was used to capture ice shedding from the 
blades. The movies were taken at 2000 frames 
per second and capture approximately a 60* arc 
of the rotor rotational field. 

Powered Force Model 

The Sikorsky PFM Is a self-contained, gen- 
eral purpose rotor test rig. The load measuring 
Systems, rotor power, and control mechanisms are 
located within the model frame. The PFM can 
accept a wide range of rotor systems and fuselage 
skins. The rig is capable of accepting a fuse- 
lage and powered tail rotor, each containing its 
own balance system. For this project only the 
main rotor and Its balance were installed. The 
attachment of power, lubrication, control and 
signal lines complete the model Installation. 

The main rotor was located 1.02 m (3.35 ft) 
above the tunnel floor, about 0.10 m (0.33 ft) 
above the centerline. The U.S. Army UH-60A 
model skins were used to enclose the balance 
and rotor hardware (Figure 3). A small rotor 
head fairing contained the hub accelerometers. 
While the fuselage skins are scaled for a 2.86 m 
(9.37 ft) rotor, tunnel dimensions dictated the 
use of a 1.83 m (6.00 ft) rotor for this test. 

The rotor head used for this project 
(Figure 5) was a Sikorsky-designed and NASA 
Langley-owned general purpose, four bladed, fully 
articulated head with coincident flap and lag 
hinaes at the 8.3 percent (76 mm or 3.00 in.) 
racial station. The hub has adjustable lead-lag 
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TABLE 1. - ROTOR 

ICING CONSORTIUM RESPONSIBILITIES FOR MODEL 
ROTOR ICING TEST 


Company 

Tasks 

NASA Lewis 

Project coordination; icing test tech- 
niques development for model rotor 
icing test 

Propul s Ion 
D1 rectorate-AVSCOM 

Icing tunnel test hardware 

Sell helicopter 

PC-based safety of flight system 

Boe 1 ng hell copters 

Test plan; analytical modeling 

McDonnel Douglas 
he 1 icopters 

Composite rotor blade design and 
construction; test support 

SI korsky A1 rcraf t 

Powered force model and associated 
hardware; data acqutstlon and reduc- 
tion; test support 

Texas A&M Uni vers 1 

ty PC-based data acquisition and 
reduction system; test support 



FIGURE 2. - OH-58 TAIL AOTOft ASSEMBLY. 
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FIGURE 5. - SIKORSKY FFH ROIOR AMD FORCE BALANCE ASSEMBLY. 


dampers and for this test It had a zero pitch- 
flap coupling (63). In order to preserve a 
consistent lead-lag damping coefficient during 
the variation In the IRT temperature, the rotor 
viscous dampers were electrically heated to main- 
tain the fluid at about 10 m C . Potentiometers 
mounted on the flap, lead-lag, and pitch axes 
provide positional Information for the rotor 
system relative to the shaft. 

The PFM main rotor was driven by two Task 
3-phase variable frequency Induction electric 
motors. As Installed, the motors can rotate at 
speeds up to 8000 rpm. A 3.5 to 1 reduction 
gear Dox operates the rotor shaft speed to a 
maximum of 2286 rpm, which equates to a rotor 
tip speed of 219 mps (718 fps). Each motor has 
a continuous rating of 60 hp. Rotor speed was 
measured using an optical encoder. Model power 
was supplied by a variable frequency motor- 
generator set provided by the Aerof 1 1 ghtdynami cs 
Directorate, U.S. Army Aviation Research and 
Technology Activity. The motor-generator set 
featured a digital closed-loop rpm feedback con- 
trol, which maintained a set rotor speed under 
the power variations that occurred due to Ice 
accretion, shedding, and control changes. 

Individual elements of the model are 
described in the following paragraphs. 

Load measurement system . - The rotor forces 
and moments (except torque) were measured using 
a Modern Machine 367-301 strain gage balance. 

The balance was glmbal-mounted to the model 
frame with a soft spring attachment, which les- 
sens the possibility of ground resonance. The 
balance was electrically heated to maintain It 
at the calibration temperature of 20 "C. Rotor 
torque was obtained by measuring the reaction 
force between the gearbox (which was mounted on 
bearings) and the gearbox frame with a Revere 
Model USPi-5-B load cell. A 102 channel slip- 
ring was usee to transfer the rotating system 


signals from blade gages, pushrods, motion hard- 
ware, and accelerometers to the fixed system. 

It also provided power to the damper heaters. 

Control system . - The swashplate control 
components consisted of rotating pushrods and 
scissors, the swashplate, three electromechani- 
cal actuators, and a stationary scissors. Con- 
trol Inputs were made at the control console 
with cyclic and collective Joys tick controllers. 
The Input signals were electronically mixed by 
the controller which then moved the model actua- 
tors. A closed-loop feedback circuit monitored 
the actuator motion. Resolved shaft axis cyclic 
and flapping first harmonic coefficients were 
obtained from potentiometers located on the 
rotorhead. These values were displayed graphi- 
cally and numerically to facilitate rotor 
trlmmi ng . 

Shaft anq 1 e control . - Shaft angle for 
this test was set by the model operator using a 
remotely driven linear actuator located on the 
test rig support frame. A pitch arm transferred 
the linear motion of the actuator to the rotat- 
ing motion of the model shaft. This resulted In 
both a translational and model pitch motion for 
the rotor. The model support shaft angle range 
was from -10 to +3 # . The angles set by the 
shaft angle control were corrected by the glmbal 
deflection angle and aerodynamic wall correc- 
tions to obtain the true total shaft angle. 

Main Rotor Blades 

Ten blades were designed and built by the 
McDonnell Douglas Helicopter Company using molds 
built by Sikorsky Aircraft. The blades for the 
1.83 m (6.00 ft) diameter rotor had NACA 0012 
airfoils with a chord of 0.124 m (4.9 In.), a 
-10* linear twist, and a taper ratio of 1. The 
blade weight had tc be kept to a minimum to stay 
within rotor head stress limits at the design 
rotor speed. The weight constraint resulted in 




a thin walled hollow spar, an unsupported trail- 
ing edge, and unprotected airfoil surfaces. 

These design restrictions led to problems during 
cold temperature operations. The blades behaved 
differently, causing variations in tracking. It 
is believed that temperature cycling and the un- 
filled trailing edge area caused the tracking 
changes. Two of the blades were Instrumented 
with strain gages in the flatwise, edgewise, and 
torsional axes to provide blade loading informa- 
tion for safety of flight. The data from the 
edgewise strain gages proved to be very useful 
In Identifying shedding events. This prompted 
the addition of root edgewise strain gages to 
the noninstrumented blades prior to the second 
entry . 

The rotor blades were marked on their upper 
surfaces with the blade number and spanwlse dec- 
ade percentage marks to provide identification 
of the blade being monitored and the relative 
location of the ice that was shed. 


INSTRUMENTATION AND DATA SYSTEM 
Test Parameters 

The test parameters measured fell Into two 
main categories - those that were of research 
Interest and those that were required solely for 
safety-of-f 1 ight. The 41 parameters in the 
first category were processed by the Dynamic 
Data Aqulsltlon System (DDA$) and saved on both 
digital tape and floppy disks. Derived parame- 
ters were also computed using the DDAS and saved 
for further analysis. The test parameters used 
for analysis Include the following: 

Main rotor balance loads (6 components) 

Main rotor torque and speed 

Tunnel temperature, and static and total 
pressure 

Liquid water content (measured by Johnson- 
Wllllams probe) 

Control positions (4) and Instrumented 
blade angles (3) 

Resolved blade flapping and coning (3) 

Instrumented blade flatwise (3), edgewise 
(4), and torsional (2) loads 

Blade root edgewise load for each blade 
(November entry) 

Glmbal and rotor head accelerations (4) 

Pushrod loads (2) 

Derived parameters calculated from the above 
test parameters will be discussed later in the 
paper . 

Each of the parameters in the list above 
(except the Johnson-HI 1 1 lams LWC and the root 
edgewise loads) were stored on an analog tape 
along with a voice track, a 1-per-revol utat Ion 
signal, and a time code. These parameters used 
40 of the available 42 tape tracks. 

Several derived parameters from the trans- 
ducer outputs were computed and displayed at a 
rate of once per second to allow the model oper- 
ator and the safety of flight (SOF) engineer to 
monitor loads and set target test conditions. 
These values accounted for tunnel blockage and 


buoyancy corrections, coordinate transform- 
ations, transfers, and balance interactions. 

The SOF system recorded an additional 47 
parameters, saved temporarily in the event of a 
model malfunction, but not archived for further 
analysis. Included on the tape were oil pres- 
sure, oil and water flow rates, drive system 
and motor temperatures, motor acce 1 erat 1 ons r 
heater voltages, limit switch and solenoid posi- 
tions, rotor speed, and reference voltages. 

Dynamic Data Acquisition System (DDAS) 

The DDAS, safety of flight system, and 
model operations consoles were all located In 
the control room of the IRT. Data parameters 
were recorded on both analog and digital sys- 
tems. The analog system was Installed primarily 
as a safety of flight system and has not yet 
been used for data analysis. The heart of the 
digital system was a Digital Equipment Corpora- 
tion POP 11/34 minicomputer. This test used 41 
of the available 128 data channels, which were 
conditioned, digitized, and transferred to the 
PDP11 by a Neff System 620. Simultaneous 
sample and hold amplifiers froze, the analog 
channels before digitizing to maintain time cor- 
relation of each parameter in each time frame 
(i.e., in each data snapshot). The system was 
configured to acquire data In both time and 
rotor domain. Data were acquired at a rate of 
16 samples per rotor revolution during the 
Icing test. Data acquisition hardware and soft- 
ware were developed to reduce the quantity of 
stored data, i.e., a continuous record of a 
150 sec icing run with a rotor tip speed of 
213.4 mps (700 fps or 37.1 Hz) would produce 
3.7 million records, more than the DDAS can han- 
dle. Therefore, data were digitized only for 
the first 10 revolutions of each second (160 
samples per data burst). The data were trans- 
ferred from the DDAS to a PC/AT computer for 
further processing. 

Safety of Flight System 

Three pieces of equipment make up the mode) 
SOF Instrumentation. An automated datalogger 
was used to monitor drive component tempera- 
tures, fluid flows and pressures, damper and 
balance temperatures, and various status volt- 
ages. Critical SOF parameters such as balance 
loads, pushrod loads, blade loads and fixed sys- 
tem vibrations were manually monitored during 
testing. Each data parameter was sent to a 42 
track FM tape system which served as both an 
Incident recorder and an analog data archive 
system. 

Data Processing Equations 

A number of commonly used rotorcraft param- 
eters can be computed from the measured test 
data. An Integral part of the processing of this 
Information is the correction of the wind tunnel 
data for the effects of the tunnel walls. These 
corrections Increase the effective tunnel speed 
and alter the rotor wake flow angles. Nondlmen- 
sional terms can also be calculated to normalize 
the data for day-to-day changes In ambient pres- 
sure, rotor speed fluctuations, and temperature 



variations. Data corrections must also be made 
for shaft torque tares, gravity tares, the 
Induced angle of attack, and three-dimensional 
buoyancy and solid blockage effects. The coor- 
dinate systems, data correction equations, and 
a description of the derived parameters will be 
Included In a future NASA Contractor Report. 

TEST PROCEDURE 


Test Techniques 

The test techniques used In the IRT were 
based on previous model operating experience and 
procedures developed during the July to October 

1988 OH-58 tall rotor entry In the the DDAS 
ence 2). Each Icing run was recorded on the DDA5 
and video systems. The video provided a viewing 
history of the Ice accretions and sheddlngs. 

There were three separate video systems: one 

for safety monitoring, one for blade tracking 
(which also provided good Ice profile shapes 
near the blade tip), and one that allowed close- 
up Images of the rotor blades. The three sys- 
tems were strobe driven by a signal off the 
rotor shaft angle encoder. This gave a ’frozen 
Image of the blade. The close-up data video 
system was Installed on a tilt and pan mechanism 
along with a 35 mm camera that had a 400 mm zoom 
lens. The data video system had the capability 
of traversing the entire diameter of the main 
rotor while allowing zoom shots of as small a 
span as 0.10 m <3.94 In.) of the blade leading 
edge. The 35 mm camera was focused on the same 
close-up viewing area as the video camera, 
allowing pictures with greater resolution and 
clarity to be taken. Both the data video and 
the 35 mm camera were triggered from the same 
strobe to provide an accurate replication o 
the video Image for the 35 mm camera. 


During the test run, Information was col- 
lected on the DDAS and the video systems, and 
the tunnel conditions were monitored. Liquid 
water content Information (rise time, cloud sta 
bllizatlon, and spray bar lag from spray Initia- 
tion) was recorded from the Control Room console 
output and a Johnson-WI 11 lams (J-W) LWC meter. 
Tables providing the main rotor operating condi- 
tions were printed out at the beginning and end 
of each run. Spray times, temperatures, and 
general comments were noted. 

Post-run information was gathered by tak- 
ing 35 mm camera photographs, recording Ice 
tracings, noting visual observations, and making 
Ice molds. Pictures were taken of the blade 
planform, an end profile, and any unusual Ice 
formations. Close-up shots of ice growths were 
taken to record their minute detail. A heated 
aluminum block with a cut-out contour of the 
airfoil shape was used to make a clean slice 
through the Ice formation. A template was then 
held against the Ice shape and a tracing Tnade. 
Measurements of the ice thickness along the pro- 
file were taken at various chord locations. 
Visual observations were recorded about the kind 
of Ice, any secondary growth, and frost forma- 
tion. Molds were taken of two blades on three 
separate occasions. 


Test Operation 

A typical test run consisted of bringing 
the PFM up to speed to exercise the rotor, then 
shutting the model down to take static balance 
and gage readings. The model was then brought 
up to operating speed and once the rotor was 
stabilized, a dynamic zero (wind off. zero col- 
lective) was taken. After this record the 
rotor rpm was dropped to some nominal level while 
the tunnel controls were set and the tunnel 
started up. When the tunnel conditions were 
stable the model operator returned to the 
desired speed, set the test point conditions, 
and took an uniced baseline point. The data 
engineer then Initiated data acquisition and 
the tunnel operator initiated the spray 
sequence. At the end of the run the tunnel rpm 
was brought down to Idle (taking approximate y 
2 min) while the model operator lowered the 
rotor speed to reduce centrifugal loads on the 
accreted ice. When the tunnel test section 
speed was below 10 kn the rotor was stopped. 
After the run the reseachers entered the test 
section and documented the results. The assem- 
bly was deiced and conditions set for the next 
run. 


TEST DATA 


Test Condi tlons 

A total of 108 test runs were completed. 

Of those, 85 were icing events with the rest 
comprised of systems check-out, gravity and hub 
tares, balancing, and baselines. The test 
matrix was set up to Include a range of condi- 
tions for a number of rotorcraft performance 
parameters under varied IRT conditions. For 
this test the temperature range was -1.7 to 
-29.7 *C, the LWC range was 0.35 to 1.24 g/rrw, 
and the MVD range was 13 to 23 pm. The predomi- 
nant condition was at -15 "C, 0.50 g/m LWC, and 
15 pm MVD. Icing encounter times were from 54 
to 168 sec. Thrust, propulsive force, advance 
ratio, and model rpm were changed, within the 
above tunnel conditions, to provide a wide scope 
nf norformance maDDina. 


Accuracy and Repeatability 

The 2o accuracy of the rotor loads was 
estimated to be: 


TABLE 2. - ROTOR BALANCE 
ACCURACY 



Balance 

Typical 


full 

measure- 


scale, 

ment , 


percent 

percent 

Lift 

0.2 

1.1 

Drag 

.2 

10 

Torque 

.7 

2 . 5 


7 




SCALING 


The rotor head motion hardware consisted of 
blade pitch, flap, and lead-lag potentiometers. 
These potentiometers were of the single turn 
conductive plastic type. The total system accu- 
racy was measured at ±0.3 # . 

The ability to repeat a test condition was 
a major goal of this test program. Rotorcraft 
icing flight test data has been subject to con- 
siderable data scatter, making the application 
of that data to code validation and basic 
research difficult. A portion of the observed 
scatter occurs due to variations in the accre- 
tion of ice on an airfoil, Reference 3, but most 
of the scatter comes from uncontrolled variation 
in the cloud. The trends from the IRT data 
examined show some scatter, but, in general, the 
data appears to be of much better quality than 
previous data available in the public domain. 
Figure 6 shows similar trends for a condition 
that was repeated four times. Figure 7 shows 
ice shape comparisons for two different repeat 
conditions and illustrates the small amount of 
variance between each of the respective runs. 



FIGURE 6. - DATA REPEATABILITY FOR FOUR REPEAT CONDITIONS 
(LWC « 0.50 g/m 3 , UVD * 15 ym. y « 0.20. OR * 206 mps, 
C L /0 * 0.06*4, T ■ -15 °C>. 



FIGURE 7. - TWO-DIMENSIONAL ICE PROFILES (IWC * 0.50 g/m 3 ' 
MVI» * 15 ym, y - 0 . 20 , OR - 183 *rs, C ./0 * 0.0GN. T ■ 
-15 D C ) . 


Similarity Laws 

In a test of this nature the concept of 
icing scaling must be addressed. Ideally, a 
complete scaling law must Include 

(1) Similarity of airflow over dry and iced 
surfaces 

(2) Similarity of blade reaction to aero- 
dynamic forces 

(3) Similarity of droplet momentum 

(4) Similarity of frozen droplet mass flux 

<5) Similarity of energy balance where 

freezing occurs 

(6) Similarity of forces which cause 
shedding 

It should be noted that this test did not 
deal with scaling in a truly rigorous way. The 
rotor chord was about 0.127 m (5 in.) to reduce 
rotor chord scale effects. The tunnel geometry 
constrained the model to a 1.83 m (6.00 ft) 
diameter rotor which results in a high (0.173) 
solidity. A full scale version of this rotor 
does not exist. A number of the above scaling 
areas were addressed and are discussed in the 
following sections. 

Airflow and Shedding Similarity 

The main airfoil parameters are the Mach 
and Reynolds numbers. The model was run at full 
scale tip speeds because blade loading and per- 
formance were important. This simulated full 
scale Mach number and advance ratio. Reynolds 
number was considered of secondary importance 
because Reynolds number effects at this scale 
and larger are generally small and an iced air- 
foil transitions to turbulent flow very early. 
The ice/airfoil Interface shear stresses (which 
dominate shedding characteristics) are matched 
at constant temperature. 

Droplet Momentum 


In order to simulate full scale droplet 
momentum It can be shown that the droplet modi- 
fied Inertia parameter must be kept constant. 
The modified inertia parameter can be expressed 
as 


yO.62,1 .62 

b 038 r 


( 1 ) 


Thus, keeping the modified inertia parameter 
constant and assuming total pressure and veloc- 
ity are the same, the ratio of model and full 
scale droplet diameters can be given as 


dm 

dF 



( 2 ) 
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If the modified inertia parameter Is held con- 
stant the local collection efficiency will be 
the same for the model and full scale (Refer- 
ence 4). Therefore, It would be desirable to 
test at very low droplet sizes, e.g., below 
10 urn. The envelope of the IRT constrains the 
droplet size to no lower than 13 pm for the con- 
ditions tested. 


Droplet Mass Flux 

It has been determined that If the accumu- 
lation parameter Is held constant, then the 
droplet mass flux which Is frozen will remain 
constant between the model and full scale. The 
accumulation parameter Is given as: 


(LHCHVHt) 

pc 


(3) 


The derivation of the accumulation parameter 
assumes : 

(1) Ice shape stays geometrically similar 

at al 1 times . , , ^ 

(2) Local collection efficiency at the 
stagnation point Is constant between model and 
full scale bodies. 

(3) Freezing fraction at stagnation point 
Is constant between model and full scale bodies. 

(4) Shedding of any water from Iced surface 
Is zero or constant for both model and full 
scale bodies. 


With the model's smaller chord, the Icing test 
time needs to be Kept relatively short In order 
to replicate a full scale Icing encounter. 


Energy Balance 

The Impinging droplets will strike the air- 
foil surface at geometrically similar locations 
as those of a full scale case If the above scal- 
ing laws are applied. However, energy balance 
concerns must be Incorporated to assure that 
these droplets freeze In the same place and man- 
ner as a full scale Ice shape. This Is accom- 
plished with the concept of freezing fraction, 
which Is the ratio of frozen liquid to the 
total amount of liquid In a control ^ vol ume and 
Includes the Impinging water as wel 1 as runback 
from neighboring control volumes. The general 
energy balance on the Interfacial layer was 
developed by Messinger and can be found In 
Reference 5. 


the similarity laws are evident. First, the 
mass and energy balance do not include the 
effects of movement and shedding of water on 
the surface as well as Ice shedding. Further, 
the basic assumptions from which the Ice accre- 
tion energy balance was derived have recently 
been shown to be in error. Films have shown 
that there Is a distinct multlzonal characteris- 
tic to airfoil ice accretion (Reference 6). The 
traditional runback model appears to be valid in 
the region near the stagnation point. Here, the 
energy balance described above would be appro- 
priate. However, In the region aft of the tran- 
sition point, which would Include the majority 
of the region of accretion, surface tension 
forces appear to dominate and little or no run- 
back occurs. Thus, more work leading to a modi- 
fied similarity law is needed in this area. 1 1 
was felt that, given the constraints of the 
project that similarity was met to the highest 
possible degree for this rotor. There lsa 
great deal of confidence that the data taken 1 s 
valid for analysis purposes. Future validation 
of this test technique however, would dictate 
that a model rotor with a full scale equivalent 
be chosen. 


DATA ANALYSIS AND CORRELATION 
Data Initialization 

The data has been reviewed to determine 
experimentally the onset of icing. The estab- 
lishment of this onset was based on a Johnson- 
Hi 1 1 lams LHC probe data when available, as well 
as the change In rotor torque and lift. This 
approach has been used throughout the analysis 
for the preparation of parameter comparisons. 

Liquid Hater Content 

The liquid water content and droplet diame- 
ter of the spray cloud are dependent on the 
spray bar water and air pressures, and the tun- 
nel velocity. The pressure values are generated 
from a set of calibration equations. The above 
approach gives an accurate setting for U*L, but 
It does not provide any history of the spray 
which Is necessary to understand the cloud time 
lag between spray on and model inundation, rise 
time to the target LHC, and cloud characteristic 
during the spray. A Johnson-Hi 1 1 1 ams <J-W) , 
uid water content meter was installed slightly 
in front of and below the rotor. This provided 
cloud information to detail the above concerns. 


Scaling theory requires that the freezing 
fraction be held constant between the scaled 
and full scale case. This requires similarity 
of the heat transfer coefficient, which is 
strongly dependent on the surface roughness. 

For this test it was assumed that, after the 
onset of icing, the surface roughness of the 
ice was similar to that of full scale case. 
Therefore, the energy similarity was approxi- 
mately satisfied. 

Due to the weight constraints the blades 
were not dynamically scaled. It should be noted 
that the current state of the art scaling laws 
are by no means all-inclusive and much work is 
needed in this area. Several shortcomings of 


The duration of the sprays ranged from 54 
to 186 sec. Scale rotor testing must use these 
shorter times to replicate full scale condi- 
tions. However, this presents some concern in 
the IRT because these short spray durations fall 
within the spray stabilization time, t.e., it 
normally takes up to a minute for the cloud to 
reach steady state condition. A plot of the 
J-H LHC data for a number of runs (Figure 8) 
indicates that there is some rise time prior to 
reaching the target value (with the rise time 
Increasing at the higher LHC runs), but the 
spray condition remains relatively close to the 
desired setting after the ramo up. The cloud 
reached the model about 10 sec after the spray 
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on" command was given. The LWC repeatability Is 
shown for two different spray bar settings In 
Figure 9. 



FIGURE 8. - VARIOUS l WC RUNS VERSUS ICING T IMF (MVD * 15 pm, 
U * 0.20. OR - 206 mps, C L • 0.064. T » -IS °C>. 
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FIGURE 9. - LWC REPEATABILITY VERSU5 ICING TIME (MVD * 15 pm, 

M - 0.20. OR • 206 mps, C|/0 - 0.06*4. T • -15 °C>. 


The data follows expected trends when the 
LWC Is increased. The torque data shows an 
Increase (Figure 10) and the lift performance 
degrades (Figure 11). 



riGiRr '0 - lw: f-fect tor a c 0 /c versus king time <kvd « 

«5 um. y -- i\?C. HR ’ 706 mps. fj /o - C.06'i. 1 - 15 °C). 



FIGURE 11. - LWC EFFECTS FOR A C^/0 VERSUS ICING TIME <MV|» - 
15 pm. y * 0.20. OR * 206 mps. C L /0 1 0.064, T « -15 °C). 


The median volume droplet diameter was 
changed in a range from 13 to 23 pm at a con- 
stant LHC of 0.75 g/nw. The results Indicate 
that an increase In MVD size has the expected 
effect of Increasing the torque (Figure 12). 

The droplet diameter has no significant effect 
on lift for the range of MVD sizes tested. 

There was good repeatability for the 13 and 
15 pm cases, but the results are less clear with 
the 18 and 23 pm sizes. 



FIGURE 12. - MVD EFFECTS FOR A Cg/0 VERSUS ICING TIME flK » 

0.50 g/m 3 . y « 0.20. I?R = 206 mps. C^/O * 0.064. T * -15 °C>. 

Delta Lift And Torque Trends 

The precision of the balance made the accu- 
racy of X bar less than desired. The balance, 
which was sized to accommodate the vibratory 
loads due to ice shedding, was too large to 
measure precisely the low drag loads at the pri- 
mary test speed of 31.1 m/sec (80 kn). However, 
the data show that shaft angle, and hence X bar, 
are of secondary importance in the data trends 
(Figure 13). Therefore, small differences in 
X bar are not significant, and data comparisons 
are made with runs that have similar but not 
identical values. 

The test was generally done with the collec- 
tive pitch and shaft angle held constant during 
the icing encounter. For most of the trends 
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FIGURE 13. - SHAFT ANGLE EFFECTS ON TOROUE (IWC - 0.50 g/m 3 , 
nvn » 15 ym. g « 0.20. OR * 206 mps. C. /O • 0.06m, T * 

-16 °C. L 



FIGURE 15. - EFFECT OF HOLDING LIFT CONSIANI TOR C L /0 VERSUS 
TIME (LWC - 0.50 g/m 3 . MVD ■ 15 pm. g - 0.20. OR - 206 mps, 
C L /0 ■ 0.06M. T - -15 °C>. 


shown In this report the Initial Cl/c was" about 
0.065 and the Initial Cq/o was about 0.0047. 
However, many theoretical comparisons may be 
better performed at constant lift conditions. 
Therefore, three conditions were repeated with 
the model pilot maintaining zero flapping and 
C{/o at the initial value. The transient nature 
or Icing creates a significant pilot workload, 
but the results show that the condition can be 
flown, and that changes in collective during the 
run have an impact on the data. Figure 14 shows 
the power Increases for this constant lift case 
and Figure 15 shows the lift history. Collec- 
tive increased about 0.5* for these cases. 

Temperature was varied from near freezing 
to very cold (-30.5 °C) . Figure 16 shows a plot 
of the power Increase for temperatures ranging 
from -3.7 to -15.2*. The torque Increases as 
temperature decreases because, in this tempera- 
ture range, the radial icing extent Increases. 

As more of the rotor blade is iced the perform- 
ance penalties become more severe. Figure 17 
shows another plot of torque rise for a much 
colder temperature range, -15.2 to -30.5 *C. 

Here, the torque rise trend reverses and actually 
decreases with temperature. As the temperature 
decreases the ice shape changes from glaze to 
rime, and the performance penalties are less. 



FIGURE *M . - EFFECT 0 r HOLDING LIFT CONSTANT FOR C Q /0 VERSUS 
T ire (LWC r 0.50 g/nv. MVD = 15 gn, y * 0.20. flP - 206 mps. 
C L /C - C.06M, T = -15 D C). 



FIGURE 16. - TEMPERATURE EFFECTS FOR A C Q /Q VERSUS ICING TIME 
(LWC • 0.50 g/m 3 , MVD « 15 gm, y * 0.20, n R - 206 mps, 
C L /0 « 0.D6M) . 



TIME, sfc 


FIGURE 17. - TEMPERATURE EFFECTS FOR ± C Q C VERSUS 
ICING Tire awe - 0.50 o/m 3 . MVP - u ' 0.2C. 

R * 206 mps. C L /0 = 0.06M-.. 
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Figure 18 shows a plot of the lift loss for the 
entire range of temperatures. The trend Is con- 
sistent with that of the torque, except for the 
extremely cold temperature (-30.5 *C). 

The effects of icing on rotor torque at 
different tip speeds is shown in Figure 19. 

There appears to be a critical velocity between 
152 and 183 mps. For the former case the torque 
rise steadily Increases throughout most of the 
run. However, the roll-off and maximum values 
for the three higher speed cases is roughly the 
same. For the higher tip speeds a balance 
between shedding and accretion appears to have 
been reached such that the torque levels off 
after roughly 50 sec of Icing time. 

Although data were taken for sweeps of 
advance ratio and disk loading at this time there 
is not enough evidence to support any conclusive 
remarks In these areas. These topics will be 
addressed in detail once a thorough analysis of 
the data 1 s completed . 



FIGURE 18. - TEMPERATURE EFFECTS FOR a C L /0 VERSUS ICING TIME 
(LVC ■ 0.50 p/m-*, MVD * 15 urn, u ■ 0.20. OR * 20G mps. 
C L /0 • 0.064). 



FIGURE 19. - TIP SPEED EFFECTS FOR ^ Cg/0 VERSUS ICING TIME 
awe - 0.50 g/m 5 . MVD * 15 ym. u e 0 . 20 . C, /0 * 0 . 064 , 

T « -15 °C>. 


Accretion 

Fundamental to any simulation of an Icing 
encounter Is the correct representation of the 
accreted ice shape. If accurate data on per- 
formance degradation is desired, a model icing 
experiment must accrete ice in a manner similar 
to full scale flight. An assessment of the 
quality of the data must begin with an evalua- 
tion of the correctness of the accreted ice 
shapes In terms of repeatability and trends. 
Repeatability has been established In an 
earlier section; a discussion of the trends 
will follow. 

A helicopter rotor blade will experience 
nonuniform ice growth as a function of radial 
location due to the velocity variation and the 
attendant aerodynamic heating effects along the 
span of the blade. The local rime Ice growth 
is dependent on the Accumulation Parameter which 
is a function of velocity, and as velocity 
increases the overall thickness of the resulting 
ice shape will increase. Also, as velocity 
increases the effects of aerodynamic heating 
become important and raise the surface tempera- 
ture to a point where rime transitions to glaze. 
These trends are born out in Figure 20 which 
shows the radial variation of accretion for a 
typical run. The primary shape increases in 
Size from the 50 percent to the 70 percent span- 
wise location, and the ice profile changes from 
a smooth rime shape at the 50 percent location 
to a rough glaze shape at the 90 percent 
location. 



CHORD, m 


FIGURE 20 - RAIMA! VARIATION OF ICE SHAPE <IWC - 0.50 C/p' ? . 
MVP = 15 ur, y - 0.2C. OR » 213 mcs, C, /C - 0.045, 
T*-25°Ci. 1 
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The accretion characteristics are also 
determined by the droplet size and the liquid 
water content. Rime Ice tends to be associated 
with smaller droplets and glaze Ice with larger 
ones. Figure 21 Illustrates this by showing the 
transition from a rime shape for 1 3 pm to a 
rougher glaze shape for 18 pm. The figure also 
shows that the lower Impingement limit moves 
back as the droplet size Increases. The smaller 
particles, with lower Inertia parameters, are 
more likely to be deflected around the airfoil 
resulting In a narrow collection zone close to 
the stagnation point (Reference 7). As the par- 
ticle size Increases the Inertia parmeter 
Increases and the droplets are less influenced 
bv the airfoil flowfleld. striking the surface 
further aft of the stagnation point. The liquid 
water content Is a measure of how much moisture 
Is In the cloud; the higher the liquid water 
content, the more severe the Icing encounter. 

As shown In Equation (3) for rime ice, the 
thickness of the Ice is linearly proportional to 
the liquid water content. This can be seen In 
Figure 22 which shows an increase in ice thick- 
ness for Increasing liquid water content. 





riGUKE 21. - ICE SHAPES FOR VARIOUS flVP'S UWC * 0.50 Q/m 3 , 
\l - 0.20. OR • 206 mps, C L /0 - 0.064, T • -15 C C r 
r/R - 501). 


The ability to adequately predict a given 
Ice shape Is critical to a rigorous methodology 
which predicts performance degradation during an 
icing encounter. The current state of the art 
In this area Is the computer analysis LEWICE, 
which was developed at NASA Lewis. LEWICE Is a 
two-dimensional code which, given the geometric 
and atmospheric conditions, will predict the 
cross-sectional ice shape for a specified icing 
time. Several comparisons have been made 
between the predictions of LEWICE and this 
experimental data. The local angle of attack 
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FIGURE 22. - ICE SHAPES TOR VARIOUS IWC'S <f1Vl> « 15 ym 
M - 0.20, OR - 206 mps, C^/O - O.Ofc'i. T * -15 °C. 
r/R • 50X) . 


was determined with the Boeing Helicopter B65 
computer performance code. The Mach number and 
angle of attack were azimuthally averaged for 
Input into LEWICE, as shown to be appropriate 
by Korkan, Dadone , and Shaw (Reference 8). The 
guidelines for time step size and surface rough- 
ness were those described by Korkan and Britton 
(Figure 9). Figure 23 shows two sample compari- 
sons for a rime condition. The overall aggree- 
ment Is good with slight overpredlctlon by 
LEWICE on the lower surface. The majority of 
comparisons made to date have been for rime con- 
ditions and the preliminary results have been 
encouraging . 

In any experiment of this type It is Impor- 
tant to establish the quality of the accretion 
data because It Is fundamental to the confidence 
In the rest of the experimental Information. 

The preliminary evidence Indicates that the ice 
accretion data Is repeatable and follows the 
correct trends. Any substantial quantitative 
statement about the quality of the data will 
have to wait until a full analysis can be 
completed . 
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FIGURE 23. - COMPARISON BETWEEN EXPERIMENTAL AND THEORETI- 
CAL ICE SHAPES awe * T.oog/ff 3 , flVD - 15 ym. P ■ 0.20. 
DR » 70S mps, C L /0 * 0.064, T - -IS °C). 


Shedding 

The Ice shedding process Is a Important 
factor In a propeller or rotorcraft Icing 
encounter; this phenomena dominates the outcome 
of the entire Icing event. The shedding proc- 
ess Is Influenced by centrifugal forces, aerody- 
namic loading, and the elastic deformations of 
the blades. Shedding can pose several potential 
problems for rotorcraft. The Ice leaving the 
blades at high tip speeds has substantial energy 
and can cause fuselage damage for rotorcraft, 
especially tilt-rotor and propeller applica- 
tions. Also, when Ice sheds In a nonuniform 
fashion an out-of-balance condition can arise 
resulting in high vibratory loads. 


large fluctuation In the strain gage reading. 

There were two major shedding events for this 
run; one at approximately 25 sec and another at 
about 38 sec. For both events blade number 4 
experiences a shed first, followed by sheds from 
the other blades. Rotor torque, also shown In 
Figure 24, exhibits a small reduction in torque 
during minor shedding, with more pronounced 
changes In torque at major shedding events. 

There were Instances where Ice was shed asymmet- 
rically, but the model vibration load limits 
were not exceeded. Figure 25 shows an example 
of a section of a blade with Ice shed. While 
this shows a relatively clean break In the ice, 
the majority of the time there Is some residual 
Ice left on the blade after a shed. 

A preliminary viewing of the high speed 
movies found nine Ice shed events either par- 
tially or completely captured. It Is hoped that 
this Information, along with the shed Impact 
data can be used to validate an Ice shedding 
computer model currently being developed. 

The condition of the blade surface has an 
effect on shedding characteristics. While the 
blades and exposed areas of the hub were wiped 
clean of any bearing grease prior to every run, 
post run examination revealed some degree of^ 
grease on the Inboard portions of the airfoil 
leading edge. The effects on adhesion proper- 
ties are not known, but the data examined show 
that Ice shedding was reasonably repeatable. 

Correlation With Theory 

The current work In analytical modeling has 
been a preliminary Investigation of ice accre- 
tion (shown In Figure 19). We will continue to 
correlate the experimental data with the B65 and 
Sikorsky GRP forward flight performance codes, 
and the LEWICE accretion code. 


Figure 24 shows a plot of the root strain 
gage output for each of the blades as a func- 
tion of time. When a shed occurs there is a 
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Concluding Remarks 

The examination of the data from the PFM 
testing has provided encouraging results. The 
quality of the test data appears to be excel- 
lent. The changes In lift and torque are well 
documented and are remarkably repeatable. The 
model Instrumentation clearly shows the time of 
shedding events and templates have been made to 
document the Ice shapes at the end of each test 
run , 

The data show the effects of temperature, 
rpm, liquid water content, and droplet diameter 
on icing. The qualitative trends of the data 
reviewed to date are correct and the test proc- 
ess looks sound. The lift at a constant control 
collective dropped by up to 15 percent. Torque 
increases of 50 percent and more were common. 

It appears that the techniques employed have 
been validated by the results obtained, and that 
the data will be useful for code and scaling 
research and development. 

Future Plans 

The near term goal is the completion of the 
data analysis for the PFM entries. We will con- 
tinue the examination of the rotor performance 
data. Only a small amount of the ice accretion 
and shedding data have been reviewed. Correla- 
tion studies will be done to see If new empir- 
ical equations need or can be developed. High 
speed 16 mm movies will be processed, and this 
Information used with the energy impact data to 
support Ice accretion and shedding model devel- 
opment. The molds taken during the testing will 
be used to make ice castings and hopefully these 
will have the quality needed for ice simulation 
experiments . 

The longer term goals Include further model 
icing tests, coordinated with full-scale rotor- 
craft icing flight testing. This would follow a 
logical progression for development and verifi- 
cation of the model rotor test techniques and 
the analytical methods, and Identify where their 
use in main rotor design is appropriate. 


It took the efforts of all the consortium 
members to make this test successful possible. 
The authors wish to thank all those Involved 
throughout the development of this project, In 
particular the NASA IRT Test Installations Divi- 
sion personnel and Sikorsky Aircraft Experimen- 
tal Aeromechanics Section model operators. 
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Appendix C 


Run No. 

LWC 

VMD 

Time 


(g/™ 3 ) 

V 

(sec) 

80 

0.98 

15 

54 

76 

0.7 

15 

200 

76 

0.7 

15 

200 

76 

0.7 

15 

200 

76 

0.7 

15 

200 

76 

0.7 

15 

200 

76 

0.7 

15 

200 

76 

0.7 

15 

200 

71 

0.43 

15 

80 

73 

0.38 

15 

80 

71 

0.43 

15 

80 

69 

0.5 

15 

168 

70 

0.5 

15 

54 

69 

0.5 

15 

168 

76 

0.7 

15 

200 

100 

1.00 

15 

80 

94 

0.5 

15 

80 

90 

0.74 

13 

168 

85 

0.74 

23 

80 

85 

0.81 

18 

80 

80 

0.98 

15 

54 

47 

0.5 

15 

80 
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